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The electrostatic potential which arises from a lattice array of point ions is computed in terms of 
a Taylor's series expansion for small distances from a lattice site. This expansion gives the change 
in electrostatic energy when an ion moves in the background of a perfect point ion lattice potential. 
The Taylor's series coefficients for terms up to fourth order in the ion displacement are evaluated for 
the NaCl and CaF 2 lattice structures. 
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1 . Introduction 

Most theoretical studies of defects in ionic solids 
require a knowledge of the manner in which the lattice 
distorts to accommodate the defect. Classical ionic 
lattice theory gives the interaction energy U^v between 
two ions fJL and v at a separation r M »>= | r M — r,,| ; namely, 

U^=(Z^Z„lr. l ,)-(C^Jrl,) 

where the four terms are respectively the coulomb 
electrostatic, dipole-dipole, dipole-quadrupole, and 
repulsive contributions to the lattice energy. The 
charge on the ion v is Z„, C^ and C ( *l are the van der 
Waals constants for ions \l and v, and the repulsive 
energy <£>rep( r /^) takes the Pauli exclusion principle 
between the /xth and the v\\\ ion cores into account. 
The repulsive energy <£ rep is a short-range function of 
r^v. The cohesive energy <J>(r ) for the crystal becomes, 



*(ro)="2 U^ + Uo, 



(2) 



where r () is the nearest neighbor distance for the per- 
fect lattice and where Uq is the lattice energy. 

The presence of a defect causes the neighboring 
ions to move from their perfect lattice sites. This 
motion will modify the energy of the lattice and this 
change in the lattice energy is most important in a 



study of the properties of defects in ionic solids. 
Among the many terms which occur in the change in 
the lattice energy due to a defect is the one which 
represents the change in electrostatic energy when 
a neighboring ion moves in the background of a 
perfect point ion lattice potential. 

As an example, let us consider the F center, which 
is an electron localized about an anion vacancy. The 
ions neighboring the F center defect move in a self 
consistent manner to accommodate the F electron. 
The many terms in the change in lattice energy are 
grouped in some convenient manner. One of these 
terms is the change in electrostatic energy AV P which 
occurs when a neighboring cation moves in the back- 
ground of a perfect point ion lattice potential [l]; 1 



AF„ 



{r l- ru)=Zv ^ Zfl {\r' v -r fl \-*-\v v -rv\-*}. (3) 



The quantity AF r (r',— r„) is the change in electro- 
static coulombic energy when ion v moves from the 
perfect lattice site r v to the position i\,, which is not 
a perfect lattice site, in the background of a perfect 
lattice point ion potential. Denoting the distortion 
(ion displacement) by r=r£ — r„ we write the energy 
of the ^th ion at the position r£=ivH-r in the form, 



F 1 ,(r)=Z„X^|r,-r M + r|-'. 



(4) 



1 Figures in brackets indicate the literature references at the end of this paper. 
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Hence, the change in the electrostatic coulomb energy 
when only the ^th ion moves in the background of a 
perfect lattice is 



A^(r) = F,(r)-F,(r=0). 



(5) 



The few researchers [2, 3, 4], who have considered 
local lattice distortions near the F center in a manner 
self consistent with the F electron state, consider 
only the alkali halides and neglect this term bV v . 
This may be reasonable in the alkali halides for 
optical absorption discussions in which the lattice 
distortions are small. But one may question their 
neglecting the term bV v for optical emission studies 
in the alkali halides and for other optical discussions 
in the alkaline earth fluorides and oxides. The lattice 
distortion may be as large as twelve percent of the 
nearest neighbor distance ro for the relaxed excited 
state in the alkali halides and in the alkaline earth 
halides and for all states in the alkaline earth oxides. 
In this paper we shall express the term AV P as a series 
expansion in the distortion r=r£ — r v and evaluate by 
Ewald's method [5] the lattice summations which 
give the series coefficients. We shall compute explicitly 
these coefficients for the NaCl and CaF 2 lattice 
structures. 

2. Formulation 

When (r/ro) < 1, we may expand the electrostatic 
energy given by eq (4) in a Taylor's series about the 
point r=0; 



V v (r) - F„(r=0) + r V V v {v) +| (r • V)*F,(r) 



+ ^(r-V) 3 ^(r)+^(r-V) 4 ^(r) + . ... (6) 

The various derivatives of V v {r) with respect to the 
Cartesian components of r are evaluated at r=0. 
Many of the derivatives in eq (6) will be zero by sym- 
metry arguments, but which ones are zero will depend 
upon the specific lattice structure and upon whether 
an anion or a cation moves. 

Let us first consider the point lattices and the re- 
ciprocal lattices for the NaCl and CaF 2 structures. 
The primitive translational vectors for these two 
lattices are 

a,= (a/2) (0,1,1), a 2 = (a/2) (1,0,1) 

and a 3 = (a/2) (1,1,0) 

and the volume of the unit cell v c is 

t>c = ai- (a 2 Xa 3 )=(a 3 /4). 

The lattice constant a is the Na — Na distance or the 



Ca — Ca distance. When a cation is at the origin, the 
position vector r/~ for the other cations (Na + or Ca ++ ) is 

r} = (a/2)[(l 2 + l :i )x+ (l, + l^y + il^ l 2 )z] 

where the /,'s are 0, ±1, ±2, ±3, etc. The anions 
for NaCl structures are located at the sites 

rr = r, + -fx„, 

where x« = (a/2) (x + y + z) and the anions for the 
CaF 2 structures are located at the sites 

rr(l) = r, + + x a (l) andrf(2) = rf + x«(2), 



x«(l)=(a/4)(i + y + z) and x„(2)= (3a/4)(i + y+£). 

We also choose the distortion r to be given by 

r= (a/4)(oii + o 2 y+o 3 i). (7) 

The position of a nearest neighbor cation when an 
anion is the reference ion is r c =(a/2)(0, 0, 1) for 
NaCl structures and r c =(a/4)(l, 1, 1) for CaF 2 
structures. The nearest neighbor distance r for the 
perfect lattice is then ro= (a/2) for NaCl structures 
and r = (V3a/4) for CaF 2 structures. Thus, 

i/ c (NaCl)=2rg and ^ c (CaF 2 ) = (16/3V3)rg. 

The first term of eq (6) is directly proportional to the 
Madelung constant for the v\h ion, a v = [V v (r= 0)ro/Z„]. 
Throughout this paper, we will use the nearest neighbor 
distance as the reference distance. 

The vectors b t = (1/a) (- 1, 1, 1), b 2 = (l/a)(l,- 1, 1) 
and b3=(l/a) (1, 1, —1) are the reciprocal triad of 
the primitive translational vectors ai, a 2 , and a 3 . The 
wave vector in the reciprocal lattice is 

g = 27r(ri]hi + rc 2 b 2 + ra 3 b 3 ) 

where the ra/'s are 0, ±1, ±2, ±3, etc. 

The series given in eq (4) is conditionally conver- 
gent. Hence, straightforward evaluations for such 
series will be most tedious and usually unsatisfactory. 
The Ewald's method [5, 6] is an elegant procedure by 
which one converts the series (4) into the sum of two 
series — each one of which converges rapidly. Referring 
the reader to references five and six for the details 
of Ewald's method, we have the following representa- 
tion for the lattice summations appearing in eq (4): 

<p(x; r)=2 |r+ + x + r|- ! 

_77 1 ^ exp (-g 2 /4G 2 ) 



te 2 /4G 2 



cos g • (x + r) 



+ ^|r + + x + r|- 1 erfc {G|r + + x + r|}, (8) 
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when 

x + r+O 

and 

^(x=0;r=0)=2|r ( + |-' 
i + o 

_77 1 v exp (~g 2 /4G 2 ) 

„. w 2- 



«eC 2 



te 2 /4G 2 ) 



2G 



+ 2|i?r 1 erfc{G|r,+ |}-— • 

I + V 77 



(9) 



We choose the quantity G so that the series in g and 
the series in 1 are both rapidly convergent. The com- 
plementary error function erfc (z) is defined by 



erfc (z) 



= (2/X^)j'~ 



exp (— y 2 )dy 



and is normalized so that erfc (0)= 1. Hence, the 
potential energy for a cation displaced a distance r 
from its perfect lattice site is 

F r (NaCl; r) = Z c [Z c <p(x = 0; r)+Z a <p(x fl ; r)], (10) 
and 
Pc(CaF 2 ; r) =Z c [Zc^(x=0; r) +Z a {?(xa(l); r) 

+ ?(xa(2);r)}]; (11) 

and the potential energy for an anion displaced a 
distance r from its perfect lattice site is 



MNaCl; r) =Z fl [Z^(x = 0; r- x„) 

+ Z«^(x«; r — x rt )], 



(12) 



and 

rfl(CaF 2 ;r)=Za[Zc^(x = 0;r-Xa(l))H-Za{^(xa(l); 



(d"V v (r)ldx*dy t dzV) for n ^ 4. The evaluation of such 
derivatives is straight forward but is very tedious. 
We tabulate in appendix A these derivatives up to 
order four; i.e., n ^ 4. 

3. Results 

Combining eqs (6), (7) and the equations in appendix 
A, we write the electrostatic energy in the form, 

VAr)=—+-U<ri + n + (ri)rn 

To To [ 



+ - {<j\ + &\ + cr :i ) Vu2 X + (o-iGT 2 + 0- 2 0r 3 + 0- :} 0-i ) V vX y 



+- (or? 4- crj + cr|)F' I ** 



+ - (erf cr 2 + or|or 3 + cr^cri + cricrf + cr 2 cr§ + ow*) Pi*** 



+ o-,o- 2 o-3^^ + — (o-J + o^ + o-^PVu- 
Z4 



+- (crio-i; + o" : JcT2 + <T2<t\ + cr:Jo- :i + aw'* + cr.^cri ) P^xy 
6 



+ - (crfcrrj + (7|(r| + <t\gt\ ) V u2 x2y 
4 



+- (erf cr 2 o-3 + <riO-!cr3 + cri <r 2 or§) F^y* + ' ' * 



(14) 



The series coefficients V v8 xtyuz are the derivatives 
{d"Vv/dx s dy t dz u ) evaluated at r=0. We list in table 1 
the Madelung constants a v and the series coefficients 
V V sxtyuz for the NaCl and CaF 2 structures. 



r— X«(l) ) H-^(Xa(2); r — X a (2))}J. (13) TABLE 1. The Madelung constant a v and the series coefficients 

V„ sxtyuz which appear in eq {14) for NaCl and CaF 2 structures 

The reference ion moves a distance r given by eq (7) from its perfect lattice site. The 
reference distance is the nearest neighbor distance r<»; r =(a/2) for NaCl structures and 
r<> = (V3a/4) for CaF 2 structures. The number /V X 10 +n is denoted by N +0n. 



The ionicity of the cations is Z c and the ionicity of 
the anions is Z (l . 

We then see from expansion (6) that the term 
F„(r = 0) and the following types of derivatives 



(d"V l ,(r)ldx°dy t dz u )\r=o 

must be evaluated, where n, s, t, and u are zero and 
positive integers and where n = s + t + u. We shall 
use the Ewald representation to compute the term 
K„(r=0) in expansion (6). Because the spatial deriva- 
tives of the lattice summation (4) for V v (r) lead to 
expressions which converge more rapidly than the 
lattice summation (4) itself converges, we find that 
Ewald's method is not necessary whenever we evaluate 
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4. Appendix A: Derivatives of the Electrostatic 
Potential 

The r • V operator which appears in eq (6) is 

_ ( , d , , d , , d\ d 4 T W5x 2 y 2 15(* 2 + y 2 ) 3 
r ' V= o-; — +cri— +oi— , (Al) 73^1=—^ 3 +> 



3* 'dy ,4 <W v ' dx 2 dy 2 r 9 r 1 r* 



where or^= (a/4)<xi. Let us denote a given term in the 
series (4) by T= (x 2 -\- y 2J r z 2 )~ li2 = r~ l . We then list 
the spatial derivatives up to order four: 

3T x 





d 3 T_ 
dx* 


15x 3 , 9% 
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